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Abstract We summarize the principles and fundamental ingredients of evolution- 
ary synthesis models, which are stellar evolution, stellar atmospheres, the IMF, star- 
formation histories, nebular emission, and also attenuation from the ISM and IGM. 
The chapter focusses in particular on issues of importance for predictions of metal- 
poor and Population III dominated galaxies. 
^ J We review recent predictions for the main physical properties and related observ- 

ables of star-forming galaxies based on up-to-date inputs. The predicted metallicity 
dependence of these quantities and their physical causes are discussed. The pre- 
dicted observables include in particular the restframe UV-to-optical domain with 
continuum emission from stars and the ionized ISM, as well as emission lines from 
H, He, and metals. 

Based on these predictions we summarize the main observational signatures (emis- 
sion line strengths, colors etc.), which can be used to distinguish "normal" stellar 
populations from very metal-poor objects or even Pop III. 

Evolutionary synthesis models provide an important and fundamental tool for stud- 
ies of galaxy formation and evolution, from the nearby Universe back to first galax- 
ies. They are used in many applications to interpret existing observations, to predict 
and guide future missions/instruments, and to allow direct comparisons between 
state-of-the-art galaxy simulations and observations. 
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Evolutionary synthesis models, first pioneered by iTinslev (Il968[ll980h . are a sim- 
ple, but fundamental tool to predict the emission of integrated stellar populations, 
such as those of distant galaxies. Their basic immediate objective is to predict the 
total spectrum or the spectral energy distribution (hereafter SED) emitted e.g. by a 
galaxy or by another ensemble of stars. 

Synthesis models are generally used to interpret observations of integrated stellar 
populations, i.e. to infer their physical properties - such as the total stellar mass, star 
formation rate, age, attenuation etc. - from comparisons between models and obser- 
vations. Often synthesis models are also used to predict/guide future observations, 
since from our knowledge of star formation, stellar evolution, and atmospheres, one 
is able to predict a large number of observables for a very broad range of parameters 
(ages, star formation histories, stellar initial mass functions, metallicities, redshifts 
etc.). This is particularly true in the present context related to the first galaxies, 
which are at the limit or beyond the reach of present-day facilities, and where pre- 
dictions are needed to plan future missions and devise observational strategies to 
search for these distant, "exotic" objects. 

Focussing on emission from stars and the ISM in the ultraviolet, optical, and 
near-IR domain (taken in their rest frame) - the spectral range where usually stars 
dominate the integrated emission of galaxies - we will describe and summarize the 
main ingredients (inputs) of synthesis models, as well as the basic assumptions and 
related uncertainties (Sect. |2]i. In Sect. |3]we show how and why physical proper- 
ties of stellar populations and corresponding observables vary with metallicity, and 
which main differences are expected between "first galaxies", metal-free (PopIII), 
metal-poor, and present-day stellar populations. Methods used to distinguish/select 
observationally PopIII-dominated and similar objects from "normal" galaxies are 
reviewed in Sect. |4] Brief conclusions are given in Sect. |5] 



2 Synthesis models - basic ingredients and assumptions 



Popular, widely used syn thesis models include the models of Bruzual and Chariot 

( Bruzual & Charlotl2003b. the Starburst 99 models specialized for young stellar pop- 

ulation s/starbursts (iLeitherer et al.119991) . the PEGASE models of Fioc & Rocca-Volmerange 
1999h. and also recent models inclu ding a special treatment of TP-AGB stars 



( Marastonl2005 ; Maraston et al.l2006l) . Recent reviews on synthesis models, includ- 



ing some basi cs as wel l as top i cs for c urrent and f uture im prov ements, have e.g. been 



presented by iBruzuall ( I2003L 1201 11) . iMarastonI ( 1201 lal) . and ILeitherer & Ekstrom 
(1201 lb . Papers presenting tests and confrontation of synthesi s model predictions 



with basic observa tions of young stellar populations include e.g. Ci d Fernandes et al. 
[2001i) : .Bruzual (.2003l) : lcid Fernandes & Gonzalez Delgado (,201oiK 



Schematically, the following needs to be known (or assumed) to be able to predict 
the spectrum of integrated stellar populations: 
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1. Stellar evolution: A description of the evolution (in time) of stars in the HR- 
diagram as a function of their initial mass, metallicity (chemical composition), 
and other parameters which may govern their evolution (e.g. initial rotation rate, 
magnetic field). 

2. Stellar atmospheres: A description of the emergent spectrum (over the spectral 
range of interest) of individual stars at all phases of their evolution. 

3. The stellar initial mass function (IMF), which determines the relative distribu- 
tion of stars of different masses at the time of formation. 

4. The star formation history (SFH) of the galaxy, describing the history of the 
amount of stars (commonly expressed in mass formed per unit time), i.e. the star 
formation rate (SFR) as a function of time. 

5. Nebular emission, i.e. the emission from Hll regions nearby massive star- 
forming regions, which - in general - cannot be separated from the stellar emis- 
sion. 

6. Attenuation within the intervening interstellar medium (ISM) of the galaxy. For 
a simple prescription this implies that we need to know the attenuation law (i.e. 
it's dependence on wavelength), and the amount of attenuation at a given refer- 
ence wavelength. 

7. Intergalactic medium opacity: Finally, since the photons emitted by a distant 
galaxy also travel through the intergalactic medium (IGM), its transmission prop- 
erties must be known/specificed. 

Items 1-3 describe the properties of stars and ensembles thereof. So-called sim- 
ple stellar populations (SSPs), corresponding to an ensemble of stars formed at the 
same time, represent the basic units. Predictions for SSPs are widely distributed in 
the literature. In practice SSPs may represent stellar populations of stellar clusters 
where the age spread between stars is small. For any arbitrary, given star formation 
history (4), the predicted spectrum can be derived from SSPs by convolution. Neb- 
ular emission (5) is important to properly describe star-forming galaxies, where the 
contribution from young massive stars is significant. The remaining items (6, 7) de- 
scribe the way the emitted spectrum is altered both at the galaxy scale and through 
the IGM on its way to the observer 

Evolutionary synthesis models traditionally describe the emission from stars 
(plus emission from surrounding Hll regions in some extensions). By construction 
such models are thus usually tailored to the (rest-frame) UV-optical-near-IR part 
of the electromagnetic spectrum, where emission from stars (H-nebulae) dominate. 
This is the domain on which the present text is focussed. Extensions of these mod- 
els to other wa velength domains, e.g . to X- r ays or to the IR-radio , have also been 



constructed by Mas-Hesse & KunthI ( 1991 ): Cervino et al. ( 2002 ). Since focussed 



on primeval galaxies and on their observability, we will here emphasize in partic- 
ular very low metallicities (necessarily more important in the early Universe) and 
relatively massive stars, which dominate the rest-frame UV emission in strongly 
star-forming galaxies. 

Let us now briefly discuss each of these "ingredients" of evolutionary synthesis 
models with a special emphasis on first/distant galaxies. 
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Fig. 1 HR-diagram for metal free (Z = 0, solid and long-dashed lines) and low metallicity (Z = 
1 / 50Z.: , dotted) stars. Isochrones of 2 and 4 Myr for Z = tracks without mass loss are also plotted 
(long-dashed). The short dashed high mass tracks evolving blueward of the ZAMS are computed 
assuming strong mass loss. The position of the ZAMS at solar metallicity (Zq) from Schaller et 
al. (1992) is sh own by the solid line. Other tracks and isochrones are also shown (see legend in 
ISchaereill2002h . Note the important shift of the ZAMS to high T ^s from low met allicity to Z = 0, 
as well as the rapid redward evolution of the massive stars. From lSchaereJ J2002h . 



2.1 Stellar evolution 



Stellar evolution models have seen a "boom" in the 1990ies, when new radiative 
opacities were published, triggering extensive evolutionary track calculations for a 
wide range of stellar masses and metallici ties. The most widely used of these mod- 
els are the Geneva and Padova tracks (see lSchaller et alJll992tlMeynet et al.ll 1994c 
Bertelli et alJll994l) . which have been extensively used since then for evolution- 
ary synthesis models and other applications. An illustration of such tracks for zero 
metallicity, and a comparison with 1/50 Z© and solar metallicity (Z©) is shown in 
Fig. [T] One of the main, if not the main property distinguishing PopIII stars from 
others is clearly apparent from this figure: the fact that massive {M>5 M^) stars 
are much hotter than their counterparts at non-zero metallicity. For the most mas- 
sive stars their effective temperature can reach up to ^ 10^ K on the zero-age main 
sequence (ZAMS). These differences has several important observational conse- 
quences, which are discussed below. A detailed discussion of the pecu liarities of in 



teri or models of PopIII stars and the literature before 2001 is given in lMarigo et al 

tooA 
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Since the stellar evolution models computed in the 1990ies including numerous 
physical processes among which in particular mass loss, large efforts have been 
undertaken to describe other processes affecting the interior evolution of stars, such 
as stellar rotation, magnetic fields, a nd the various transport mechan isms related to 
it. For reviews on these issues see e.g. iMaeder & Mevnell ( l2000l l 1201 21) The impact of 
other, maybe more "exotic" phenomena, like possible variations of the fine-structure 
constant on stellar p roperties at very low metallicity have e.g. been explored by 
Ekstrom etaD(l2010h . 

Several arguments, both observational and theoretical, indi cate that the impor- 
tance of rotation becom es stronger for low metallicities (see e.g. Maeder et al.ll 19991 : 
Chiappini et al.l 120061) . However, the predicted stellar tracks depend strongly on 



the initial rotational velocity (a free parameter) and its history (which is predicted 
from the model following also the evolution of angular momentum), which are 
very uncertain, and difficult to constrain at (very) low metallicity. In consequence, 
the impact of stellar rotation on predictions for integrated stellar populations re- 
mains currently poorly known. Preliminary explorations of t he impact of st e llar ro - 
tation, mostly at solar meta lhcity, have been presented by IVazquez et al.l (l2007h . 
Leitherer & EkstromI (1201 ih . and Levesque et al. (2012, submitted). 

Recent PopII I tracks including the effec ts of stellar r otatio n and magnetic fields 
are presented in Ekstrom et al.l (l2008l) and lYoon et al.l (|2012|) . For the initial rota- 
tional velocities chosen in the former paper, rotation has a small impact on the evo- 
lution in the HR-diagram and hence on predictions from synthesis models. How- 
ever, the surface abundances of these stars are significantly modified, with several 
implications on nucleosynthesis and chemical evolution. Effects with potential im- 
pact on evolution ary synthesis models can be found for very high rotation rates (cf. 
Yoon et al.ll2012l) . when stars are very strongly mixed, following a nearly homoge- 
neous evolution, which implies much hotter temperatures and a blue-ward evolution. 
In this case properties such as their UV and ionizing flux are significantly altered, 
leaving imprints on the spectra predicted from stellar populations containing such 
stars. However, the distribution of rotational velocities remains poorly know, es- 
pecially at very low metallicity, and hence the proportion of sta rs with properties 
significantly altered by rotation are not known. ISchaererl (l2003b has explored in a 
simple way the possible impact of very hot, homogeneous PopIII stars on the hard- 
ness of the ionizing spectrum of stellar populations. It is quite evident that we cur- 
rently do not have a clear view of the impact stellar rotation and magnetic fields may 
have on evolutionary synthesis models, quite independently of metallicity. Progress 
in these areas is ongoing. 

Other issues of general importance for synthesis models include e.g. the impor- 
tance of thermally-pulsating asymptotic giant branch (TP-AGB) stars, whose con- 
tribution to the integrated light is being debated (see Maraston 2005 . 201 Ibl and 
references therein), and the recurrent question of the importance of binary stars (cf. 



Vanbeveren et al.ll998l:ISchaerer & Vaccall998l:lBelkus et al.boolloionne & Robert 



2006HEldridgell201 1|) For our present objective, first galaxies, uncertainties related 
to TP-AGB stars are not a major concern, since we are dealing with young stellar 
populations and focussing on the UV to optical (rest-frame) domain, where these 
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Stars can safely be neglected. The evolution of close massive binary stars can al- 
ter the predicte d spectra of young stellar populations, as illustrated e.g. recently by 
Eldridgd (1201 lb . However, the physics of these stars is even more complex than that 
of single stars, and depends on a number of additional, poorly constrained parame- 
ters. In any case, it is important that the predictions from stellar evolution models be 
compared to and tested against direct stellar observations in our Galaxy and in the 
nearby Universe. This includes e.g. observations of individual stars in clusters, their 
basic stellar parameters and surface abundances, color-magnitude diagrams, inte- 
grated colors of clusters, statistics of stars of different types and trends with metallic- 
ity etc. Such compariso ns have e.g. been carried out with the (non-rotating) Padova 
and Geneva models (cf. Maeder & Mevnet|[l99i: lBruzualll200l iBruzual & Chariot! 
20031: lMassevll2003l:lMevnet & Maedeill2005b . and are being carried out for the lat- 
est tracks including stellar rotation. They serve thus also to "calibrate" uncertainties 
or unknowns in the stellar models, and hence to place synthesis model predictions 
on the best ground. 



Fig. 2 Ionising photon flux 
per unit surface (top panel) 
and liardness of tlie ionising 
spectrum (middle and bottom 
panel) as a function of ef- 
fective temperature predicted 
from non-LTE plane parallel 
TLUSTY Pop III models of 
various T^ff and logg (solid 
black lines). Triangles show 
calculations from the spher- 
ically expanding CMFGEN 
models for Pop III stars. 
Green dashed lines show 
predictions from blackbody 
spectra. Similar to Fig. 2 from 
ISchaereJ f2Q03) . 







25 
20 

m 

2 15 






^X,^\ q(He) ; 


10 


'■ — H 


H \ 1 1 1 1 1 


q(He')\ J 
1 1 \ \ — '■ 









, ^He)/q(H) J 


IT 












q(Ile*)/q(II) V 




o 

-4 






1 — 1 — — I — \ — 1 — 1 — — : 


1 

0,8 


; H 


H \ 1 \ \ \ \ 




0,6 






q{He)/q(H) J 


tf 0,4 




q(He*)/q(H) 




0,2 













i~-J^ 


< ll 1 \ 1 L_ 




150 100 

Tetf [kK; 


50 



2.2 Stellar atmospheres 

It is well known that strong departures from local thermodynamic equilibrium (LTE) 
occur in the atmospheres of hot stars. To properly predict the emergent spectra of 
massive stars it is therefore essential to use non-LTE model atmospheres. This state- 
ment holds in particular also for low- and zero metallicities, as departures from LTE 
significantly alter the level populations of H and He, which - together with electron 



Evolutionary synthesis models as a tool and guide towards the first galaxies 



7 



scattering - are the main opacity sources in stars of such composition. While gener- 



ally s tellar winds also affect the predicted spectra of massive star fe.g. lGabler et al 



1989I) . the mass outflow from very low metaUicity stars is very low, and its effect 
can hence be safely negl ected to compute the observable properties of PopIII an d 
similar massive stars (cf. ISchaereiil2002t lKudritzkill2002l : iKrticka & KubatI 120061) . 
Using non-LTE plane parallel model atmospheres is therefore sufficient for this pur- 
pose. For less extreme metallicities, appropriate model atmospheres describing main 
sequence stars , Wolf-Rayet star s with strong winds, and also cooler stars must be 
combined (cf. Conti et al] 20081) t o ach ieve the most reliable synthesis model pre- 



dictions (e.g. lLeitherer et al.lll999l:lSmith et al...200Z) 



Fig. 3 Comparison of input 
stellar spectra for photoion- 
ization models showing a re- 
alistic SED (blue line; Pop III, 
Salpeter IMF 1-100 M,;, at 
zero age) and a black body 
spectrum (red line) with the 
same average photon energy 
(£ = 30.033 eV) in the Lyman 
continuum. The SED model 
is normalized at 1240 A and 
rescaled to match the same 
ionizing flux Q{H). Despite 
having the same Lyman con- 
tinuum flux and average mean 
energy, the black body differs 
significantly, e.g. in the He+ 
continuum (A < 228 A), and 
in th e observable UV. Fig ure 
from lRaiter etai] ilOWm . 
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The predicted ionizing fluxes of H, He, and He+ from PopIII stars as a func- 
tion of effective temperature is shown in Fig. |2] For comparison, the same values 
predicted using simple black-body spectra are also shown. While for Tgff > 40 kK 
black bodies provide a good approximation to the total number of H ionising pho- 
tons (i.e. at energies > 13.6 eV), this is not the case for the shape of the spectra 
at higher energies, for the number o f He and He^ i onizing photons, and for other 
features such as the Lyman break (cf. Schaererll2002 ). Also, using black-body spec- 
tra is a bad approximation to compute observable properties such as UV fluxes and 
emission line strengths (equivalent widths), since black-bodies do not reproduce the 
main bound-free edges, which significantly shape true stellar spectra. See Fig. [3] 
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for an illustration. In short, appropriate state-of-the-art stellar atmosphere models 
should be used for rehable predictions both of individual stars and integrated stellar 
populations. 



2.3 IMF 



As for other ingredients summarized above, the stellar initial mass function (IMF) 
is a separate topic, to which many papers and conferences ar e devoted (see e.g. 
Gilmore&Howellll 19981: iBastian et al.feoKA iTrever etaLlboilb 

In the nearby Universe the IMF is well described b y a log-normal function for 
masses below ~ 1 Mq, and by a pow er-law above tha t dChabrieilbOOll) . as already 
found e.g. in the pioneering stud y of Salpetei (Il955h . Other authors approximate 
the IMF by piecewise power laws ( lKroupal200lb . In our Galaxy and in oth er nearby 
systems it is found that the IMF is only weakly dependent on environnent ( IChabriei 



20031; iBastian et al.ll2010l) . For obvious reasons it is not possible to constrain the 



IMF in the same manner (e.g. by direct star counts) outside the local Group. Hence 
our empirical knowledge on the IMF and possible variations with metallicity, stellar 
density, UV radiation field and other factor which may influence it, is very limited. 

For many/most observables measurable for high-redshift, star-forming galaxies 
only the upper part of the IMF is relevant, since massive stars largely dominate the 
luminosity in integrated stellar populations. As a reasonable rule of thumb one may 
therefore simply rescale the results derived for one given IMF to those expected for 
another IMF differing in the domain of low stellar masses (i.e. typically at M ^ 5 
Mq). This is e.g. used to "correct" results from synthesis models computed for a 
Salpeter IMF extrapolated down to 0.1 M^^, to those expected for a more realistic 
IMF below < 1 M,:;. However, the applicability of this "rule" depends on the age 
and star-formation history, and on the observable used. 

Several studies, following different arguments have suggested deviations from 
the "local" IMF in distant galaxies with deviations depending on environment, 
galaxy type, gas density, metallicity and other factors. For an overview over this 
vast topic the read er is referred to a recent conference proceedings on the IMF 



(ITreyer et al .11201 Ih . For metal-poor environments and Population III stars, simu- 



lations have long suggested a preferential mass scale much higher than in the local 
Universe, with "typical" masses of the order of ^ 10-100 Mg,. See the chapters 
of Glover and Johnson in this book, and also recent conference proceedings on the 
first stars (e.g. Hunt et al. 2008 : Whalen et al. 2010l) . This difference, basically due 
to significantly reduced cooling of metal-poor or fr ee gas, is predicted to occur be- 
low a certain critical metallicity Z^nt < 10^^^' Zq dSchneider et al ] |2002ll2003h . If 
correct, this implies a different, more massive IMF, at least below this metallicity 
threshold. Accor dingly variou s parametrisations of the IMF have be en used in this 
domain (see e.g. iLarsonll 19981: iTumlinson 120061: iRaiter et al.ll2010bh . However, as 
already mentioned above, our knowledge of the IMF, especially in such extreme 
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conditions remains very limited, and synthesis models can simply assume different 
cases and examine their implications. 

Synthesis models also generally make the assumption of a continuous, well- 
populated IMF, which is correct for stellar systems with a large enough number 
of stars. However, in "small" stellar populations the analytic statistical descrip- 
tion of the IMF will break down, and sampling of the IMF with a small/finite 
number of stars may lead to significant differences, due to the progressive ab- 
sence of massive stars. Examples of such "stochastic" IMF effects on colors, ion 



izing fluxes and others h ave been illustrated by various authors (see ICervino et al 



2000t Cervino & Luridi ana 2006; Bruzual 2002, and references therein). Typically 
stochastic effects are found to be significant for cluster masses below < 10^ Mq or 
for low star formation rate s SFR < 1 M r: yr^^ although this li mit ag ain depends on 
the observable of interest (ICervino et a l. 2000; Fumagalli et al. 201 lb . For small en- 
tities, such as the first galaxies, stochasticity may thus be important for some cases. 

The question of a stochastically sampled IMF is also related to the concept of the 
integrated galact ic initial mass function (IGIMF') introduced by Kroupa an d collab- 
orators (see e.g. IPflamm-Altenburg et afll2007t IWeidner & Kroupall2006 ). Various 
ap plications of synthe s is models implementing a st ochastic IMF can e.g. be found 
in (Pindao et all (l2002h . Ipumagafli etall (l201ll) . and lEldridgd (l201ll) . For example, 
a stochastic IMF leads to a significant scatter in the relative Ha and UV output for 
populations with a low SFR (i.e. forming a small number of stars), as illustrated 
in Fig. |4] This is due to the fact that these two emissions originate from stars of 
somewhat different mass regimes. Ha being due to more massive stars than the 
UV continuum. Since stochastic effects are of increasing importance for low mass / 
low SFR galaxies, this issue should be relevant for studies of the first galaxies. The 
impUcations remain largely to be worked out. 



2.4 Star formation history 



Star formation histories are a key ingredient for evolutionary synthesis models. 
As already mentioned above, simple stellar populations (SSPs) represent the ba- 
sic units, and for any arbitrary star formation history (SFH), the integrated spectrum 
can be derived from SSPs by convolution. Historically, simple parametrisations of 
the SFH have been used, the most common one being a family of exponentially de- 
creasing SFHs with SFR(r) <^ exp(— r/r), where t > is a characteristic timescale. 
It is well known that such star formation histories are able of reproducing the ob- 
served spectro-photome t ric properties of present-day g alaxies of all Hubble types 
( Larson & Tinslevl[l978 ; Tinslevl l 1968HKennicut3ll998h . 

Constant SFR, corresponding to T = oo, is a limiting case often used to derive cal- 
ibrations for the SFR from various observable s (e.g. UV c ontinuum, H recombin a- 
tion lines, bolometric luminosity etc.). See e.g. lKennicuttl(ll998h ; IS chaerer (l2000l) . 

Depending on the application, other, more complex SFHs have been explored. 
Many numerical simulations of galaxy formation and evolution (semi-analytical 
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Fig. 4 The ratio of SFR 
measured by Ha and UV 
fluxes versus the SFR from 
UV flux. The aster isks are 
the obs ervations of iLee et"ai] 
j2009j) , while the shaded 
region show the density of 
our individual realisations 
of synthetic galaxies. The 
thick solid lines indicate the 
mean ratios for the synthetic 
galaxies and their Iff limits. 
The dashed lines show the 
mean ratios for the other IMF 
filling method with the same 
stellar population. The upper 
and lower panels distinguish 
a stochastic sampling of 
the IMF and a sampling 
including a description of 
the cluster mass function. 
Left/right panels distinguish 
synthesis model including 
single/(single+binary) stars. 
From iEldridge aol% . 
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models, hydrodynamic models and others) are now "coupled" with evolutionary 
synthesis models to predict observable properties consistent with their (complex) 
star formation histories. 

For distant (z > 2) galaxies SF histories and corresponding timescales are cur- 
rently vividly debated. For example, simulations suggest rapid growth of galaxies 
with increasing star formation rates during the first Gyr of the Universe, i.e. at z > 6 
(^Jinlatoret al. 20071 120T Ij). From the apparent tightness of the mass-SFR relation 
of galaxies at high redshift s ome authors argue for rising star formatio n histories 
for galaxies down to z 2-3 ( Maraston et al.ll20"lot Papovich et al.ll2011 ). At some- 
what lower redshift, other authors suggest e.g. "delayed" star formation histories 
SFR(f ) 0= t exp(— f /t) showing both phases of increasing (for f < t) an d decreasing 



SFR i t > t), and suggest relatively long timescales of several 100 Myr (iWuvts et al 
From the clustering of z 4-5 Lyman break galaxies (LBGs), iLee et al 



(l2009bh argue for relatively short duty cycles (< 0.4 Gyr). Even shorter timescales 
and declining SFHs are favored from SED analysis of Lyman break gal axies at z ^ 
3-6 w ith spectral templates accounting for nebular emission (Schaerer & de Barros 
201 ll de Barros et al. 2012, in preparation). The question of star formation his- 



tories of distant galaxies is closely related to the debate about the main mode of 
star formation, i.e. schematically about the relative importance of (cold) accretion 
driven star formation and mergers (see e.g. iKhochfar & SiUduOl ll and references 
therein) In any case the question of the "typical" star formation history of distant 
galaxies is not yet settled, and a large diversity of histories - more complicated than 
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simple, parametrised functions - must occur in nature. Furthermore it is clear that 
predictions from evolutionary synthesis models strongly depend on this important 
quantity. 



2.5 Nebular emission 

Stars (and stellar populations) described above provide a dominant, but not the 
sole source of emission in the (rest-frame) UV, optical, and near-IR light of galax- 
ies. Emission from ionized regions (the so-called Hll regions) of the surrounding 
ISM is another important contribution, which needs to be taken into account to de- 
scribe/predict the spectra or SEDs of galaxies in this part of the electromagnetic 
spectrum. Hence this should, in general, be treated (added) to the prediction from 
"standard" evolutionary synthesis models describing stellar emission only. 

Most galaxy types across the Hubble sequence show at least some signs of emis- 
sion lines in the optical. Their strength increases towards late types and irregular 
galaxies, b asically due to an increasing ratio of present over average past star for- 



mation (cf. lKennicuttl|1998} . Nearly by definition, nebular emission is associated 



with star-forming galaxies, since - as long as stars more massive than 5-10 M(; are 
formed - star formation always implies the emission of UV photons with energies 
> 13.6 eV capable of ionizing H and other elements in the ISM, causing a plethora 
of recombination lines. In addition several processes in Hll regions, including re- 
combination and 2-photon emission from the 2 ^S level of hydrogen, produce a 
continuum emission longwar d of Lya whose emissivity increases with wavelength 



( Osterbrock & Ferland 20061) . Hence, nebular emission implies a priori both line 



and continuum emission. An illustration of the three emission components, stars, 
nebular lines, and nebular continuum, observed in metal-poor, nearby galaxies is 
shown in Fig.|5] 

To first order, the strength/luminosity of both lines and nebular continuum emis- 
sion depends on the flux/luminosity of ionizing photons in the Lyman continuum 
denoted here by 2(H), as can be derived from well known nebular physics, as- 
suming e.g. so-called case A or B and typical nebular densities and temperatures 
dOsterbrock & Ferlandll2006l) . n this approximation it is straightforward to predict 



the emission from the major H and He recombination lines and continuum emission 
for a given source spectrum, such as calculated by evolutionary synthesis models. 
For other emission lines from Hll regions, mostly forbidden metal-lines (in the opti- 
cal domain), full photoionization models need to be computed, or other prescriptions 
be used. 

To second order nebular emission depends on the conditions in the nebula/ISM, 
which are primarily described by the electron temperature and density (Tg, vie), and 
most importantly on the ionization parameter U . The latter depends on the radiation 
field, geometry, and density, and in general the nebular conditions also depend on 
metallicity. It is well known that Hll regions and galaxies at low metallicity show 
higher excitation and stronger optical emission lines, as illustrated e.g. in nebular 
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Fig. 5 Observed spectra of 
two metal-poor Hll galax- 
ies from the sampl e of 
lOuseva et al. focussing 
on the region between the UV 
and optical domain. Super- 
posed are the stellar and 
gaseous (nebular continuum) 
contribution. The flux is given 
in Fi units between 3200 and 
5200 A with spacing of 100 
A between tickmarks. The 
Balmer jump (not a "break" 
here), due to nebular free- 
bound emission, is clearly 
seen around 3650 A; no stel- 
lar source is known to show 
such emission. Figure ada pted 
from lGuseva et all )2006h . 




diagnostic diagrams ( Baldwin et al.ll981 : Veilleux & Osterbrockll987 ). an example 
of which is shown in Fig.|6] 



Fig. 6 The [Nil] /Ha ver- 
sus [OlII] /Hj8 diagnostic 
diagram. Small grey points 
show ~ 96000 objects from 
the Sloan Digital Sky Survey, 
the large black circle the pe- 
culiar, low-metallicty z ~ 2 
gal axy BX41 8 disc overed 
by lErb et alj i l2010l) . Col- 
ored points are predictions 
from photoionization models, 
coded by ionization parameter 
as labeled at low er left. From 
lErbetal.l i [20lo[) . 
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Since both the excitation of the ionized gas and the ionizing flux increase 
with decreasing metallicity, nebular emission must on average be stronger in un- 
evolved/early galaxies than in present-day galaxies. How this affects various ob- 
servables will be discussed in detail below. 



Indeed, for Population III stars and for ensembles thereof, ISchaeren ( 120021) has 
demonstrated that nebular continuum emission dominates the spectrum at wave- 
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Fig. 7 Spectral energy dis- 
tribution of a young (zero 
age) PopIII starburst with 
an IMF extending to 500 
M,v:. The SED shown in full 
lines includes H and He re- 
combination lines, nebular, 
and stellar continuum. The 
pure stellar continuum (ne- 
glecting nebular emission) 
is shown by the dashed line. 
For comparison the SED of 
the Z = 1/50Z;., population 
is shown by the dotted line. 
The vertical dashed lines indi- 
cate the ionisation potentials 
of H, He", and He+. Note 
the presence of the unique 
He II features (shown as thick 
dashed lines) and the impor- 
tance of nebular continuo us 
emissi on. From [Schaererl 
<2002). 



Pop 111: Salpeler IMF (1-500 M^) 




2000 4000 6000 
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length longward of Lya, as illustrated in Fig. |7] At higher metallicities the effect 
can also be significant, as discussed below. In some metal-poor, nearby galaxies, 
nebular continuum emission is readily observable in the UV-optical domain (see 
Fig. |5]l, and signs of nebular emission are frequently observed in spectra and pho- 
tometry of starburst galaxies, e.g. in diffuse, outer regions, or in the near-IR do- 
main in young star-forming galaxi es such as Hll galaxies, blue compact dwarfs 
and similar objects (cf.|lz otov et alii 1997 : Papaderos et al ]|2003l: llzotov et al ] l200lf 



Vanzi et alj|2002t iGuseva et alj 120041) . At high redshift a peculiar, lensed galaxy 
(the so-called Lynx arc at z = 3.357) has been found , whose spectrum ap pears to be 
dominated by nebular continuum emission and lines ( Fosburv et al j2003l) . A similar 
explanation has been suggested for a z = 5.5 galaxy ( Raiter et alj 2010al) . although 
alternative explanatio ns based on more conventional galaxy SEDs can also be found 
dVanzella et alj|2010[) . n any case, the objects and conditions mentioned here rep- 
resent probably just somewhat extreme examples illustrating the potential effects of 
nebular emission, and their detailed importance needs to analyzed on a case-by-case 
basis. However, there is growing evidence that nebular emission plays a role for the 
interpretati on of observations of many/ most Lyman break and Lya emitter galaxies 
at z>3 (cf. Schaerer & de Barrosll201 1 , and discussion below). 

Several evolutionary synthesis models, such as PEGASE, GALEV and others 
have long included nebular emission (continuu m emission and also lines in some 
cases, and relying on different assumptions) (cf. Fioc & Rocca-Volmerangdll999t 



Leitherer et al.l 1 19991 [Charlot & Longhettil I2OOII: lAnders & Fritze-y. AlvenslebenI 



2003 



Zackrisson et al.l2008l) . However, other widely used codes (e.g. lBruzual & Chariot 



20031) neglect this component. To interpret emission line observations of Hll re- 
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gions and galaxies, the stellar spectra predicted by ev olutionary synthesis codes 
have often been used as inpu t forphotoionisation models ( Garcia- Vargas et aljl995 



Stasiriska & Leithered Il996t IStasiriska et all I2OOII: iKewlev et alj l200lh Recentlv. 



nebular emission has again received considerable attention, in particular for the in- 
terpretation of photometric observations of dis tant Lyman break galaxies (LEGS' ) 
and Lyg emitters (LAE) as discussed below ( Schaerer & de BarrosI I2OO' 
Ono et alJlioiotlSchaerer & d e Barros 201 ll lAcquaviva et al.ll201 lb . 



201C; 



Now that all emission sources have been "assembled", we need to briefly discuss 
absorption/attenuation professes occurring along the lines-of-sight to the observer. 



2.6 Attenuation law 



To compare predictions of evolutionary synthesis model with observations, the ef- 
fect of interstellar reddening (extinction and/or attenuation) obviously needs to be 
taken into account. Simple prescriptions are generally used, describing e.g. t he mean 
obser ved attenuation law of star-forming galaxies (e.g. the "Calzetti law" ICalzetti 
200 ll) o r various laws describi ng extinction in the Galaxy or the SMC. At very high 
redshift Gallerani et al. I (l2010l) have recently argued for an attenuation law differing 
somewhat from the Calzetti law, probably due to different dust composition from 
supernovae dominating in the early Universe. 

In the nearby U niverse, s tellar emission is generally found to be less attenuated 
than emission lines ICalzettil (1200 ih . which is attributed to geometrical effects. If this 
also applies to distant star-for ming galaxies remains, however, controvers ial and 
may vary from case to case (e.g. Hainline et al.ll2009 : Yoshikawa et alj|2010l) . Other 
complication arise e.g. in the presence of multiple stellar populations suffering from 
different attenuations. More complex prescriptions and additional assumptions are 
needed to describe such cases (see e.g. lCharlot & Faill2000t l iGranato et al ] l2000l) . 

In concrete applications of synthesis models to fitting (restframe) UV-optical 
photometric observations of star-forming galaxies attenuation/extinction is often de- 
generate with the age of the stellar population, since both lead to a redder SED on 
average. This well-known age-reddening d egeneracy is ill ustrated e. g. in s tudies of 
simple stellar populat ions (see Icid Ferna ndes & Gonzalez Delgad^ 201(i). distant 
star-fo rming galaxies ( Pello et al. 19991) . distant red galaxies IPozzetti & Mannuccil 
(l2000l) . and in many other works. Obviously, the use of other, more sensitive 
indicators of age a nd/ or atte nuation, e.g. spectral lines, the Balmer break etc. 
Cid Fernandes et al.l (cf. l200lb . can reduce such degeneracies. A more detailed dis- 
cussion of this and other issues related to the fitting of observed SEDs and spectra 
is beyond the scope of this Chapter. 
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2.7 IGM 



Finally, the collective effect of the Lyman forest on predicted spectra of distant 
galaxies must also to taken into account. Indeed for redshifts above z> 3-4 the 
IGM significantly reduces the flux shortward of Lya in the restframe of the galaxy 
causing the Lyman (continuum) break at 912 A to shift rapidly to Lya 1216 A). 
Generally the average IGM attenuation is described by a s imple expressi on derived 
from statistical analy sis of the Lya forest, e.g. fo llowing iMadaul d 1 995b or recent 
updates of this work dFaucher-Giguere et al ] |2008h . 

A related, important question in several contexts concerns the transmission of 
Lya photons emitted by high redshift galaxies. Indeed, as strong Lya emission is 
one of the expected signatures of primordial and very metal-poor stellar popula- 
tions (cf. below), knowing its escape fraction out of galaxies and the subsequent 
transmission through the surrounding IGM are of fundamental importance. By the 
same token, observations of Lya emission (e.g. the Lya luminosity function of 
galaxies) and its evolution with redshift are used to constrain the IGM transmis- 
sion and hence the reionization history of the Universe. For more information on 

I — : 1 

these related topics we r efer the reader t o a vast literature including e.g . Haiman] 
(l2002h:lFan et alj (l2006h : rSchaereil(l2007l) : lHaves et all tOl lb : lDaval et aU (i201 iT : 
Lidman et al.l (120 1 2h and references therein. 



ZAMS, Z = 0-0-02. Salpeter IMF 1-100 M^, 



Fig. 8 Predicted SEDs in- 
cluding Lya and He II X 1640 
emission lines for zero age 
main sequence (ZAMS) mod- 
els at different metallicities. 
The metallicities Z = 0. (Pop 
III), 10-', lO-\ 0.0004, and 
0.02 (solar) are from top to 
bottom in the EUV (A < 
912 A), and reversed at longer 
wavelengths. The dashed lines 
are the pure stellar emission, 
the solid lines show the total 
(stellar -|- nebular) emission. 
A Salpeter IMF from 1-100 
M,.;, is assumed he re for all 
metall icities. From lSchaereJ 
j2003l) . 
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Fig. 9 Predicted model 
spectra of a PopIII galaxy 
(top panel) and a mod- 
erate sub-solar metallicity 
(Z = 1/50Zb,, bottom panel), 
at ages of 10, 100, and 500 
Myr (from bottom to top in 
each panel), evolving with 
a constant SFR. Solid and 
dotted lines show cases com- 
puted for 50% and no escape 
of Lyman c ontinuum pho tons. 
Figure from llnoudfloTlh . 
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3 From present-day metallicities back to the first galaxies 



The Starburst99 evolutiona ry synthesis models (iLeitherer et al ]ll999h and the mod- 
els of lSchaereii ( I2002ll2003h allow one to examine in detail the dependence of the ex- 
pected observational properties of star-forming galaxies with metallicities from cur- 
rent (solar) m etal Ucities to metal-poo r and metal-free cases. Furthermore ISchaerer 
(120021 l2003h and lRaiter et alJ (l2010al) also present calculations for a variety of dif- 
ferent IMFs, such that IMF or combined IMF/metallicity changes can be examined. 
We here summarize the behavior of some of the main observables based on these 
models. 



3.1 UV-optical: stellar and nebular continuum emission 

Figure|8]shows variations of the UV spectrum, including the ionizing « 912 A) and 
non-ionizing part (A > 912 A) of a young stellar population with metallicity. Three 
main features are immediately clear: First, due to the increase of the average stellar 
temperature, the stellar spectrum becomes harder/bluer with decreasing metallicity 
Z. Second, the contribution of nebular continuum emission increases concomitantly. 
At young ages, the continuum longward of Lya (> 1216 A) will be dominated by 
nebular emission. Finally, H and He emission lines become also stronger with de- 
creasing Z (cf. below). The models in Fig.|8]assume so-called case B recombination, 
i.e. in particular that all ionizing photons are absorbed within the surrounding Hll 
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region. This case obviously maximizes the nebular emission, whereas in a more 
general case where a fraction /esc of Lyman co ntinuum photons escape, nebular 
emission will be decreased; see e.g. llnoud (1201 lb and Fig.|9] 



ZAMS values 



Fig. 10 Predicted UV slope 
i3i500 for different IMFs 
(color codes) plotted as a 
function metallicity. The 
values are shown for very 
young (ZAMS) populations, 
which correspond to the bluest 
possible slopes (i.e. minimal 
j8 values). The upper set of 
lines shows the UV slopes of 
the total spectrum (stellar + 
nebular continuum), the lower 
lines using the pure stellar 
spectrum. Black lines show 
the predictions for a Salpeter 
IMF slo pe. All color codes a re 
given in lRaiter etal] jlOlObI ). 
from where the figure is taken. 
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Fig. 11 Temporal evolu- 
tion of the UV slope j8i5oo 
derived from synthesis mod- 
els of different metallicities 
and for instantaneous bursts 
(solid lines) and constant SF 
(long dashed lines). Black 
lines show solar metallicity 
models, red lines show metal- 
licities between Z = 10^' and 
zero (PopIII) and blue lines 
show intermediate cases of 
Z = 0.004 and 0.0004. The 
dotted lines show /3 if nebular 
continuous emission is ne- 
glected, that is, assuming pure 
stellar emission. Note espe- 
cially the strong degeneracies 
of j8 in age and metallicity 
for bursts, the insensitivity of 
j8 on Z for constant SF, and 
the rather red slope for young 
very metal-poor bursts. From 
[Schaerer & Pelld C2005.) . 
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The color of the UV continuum and its variations can be quantified e.g. by the 
exponent j3 of a power law adjusted to a specific spectral region. This so-called 
"j3 -slope" determined from models and measured by spectroscopy or photometry is 
often used in the literature, in particular to estimate the amount of dust attenuation 
in high-z star-forming galaxies (see e.g. iMeurer et al. I I1999I: iBouwens et al. I l2009[ 
and also the Chapter by Dunlop in this book). 

Figure [To] shows the dependence of the predicted j3 -slope (here measured be- 
tween 1300 and 1800 A) of a young (zero age) population on metallicity and on 
the IMF (coded by different colors). Fig.[TT]shows the same quantity as function of 
age for instantaneous bursts and constant SFR. From Fig.[TO]we see how the stel- 
lar spectrum steepens (becomes bluer) with decreasing Z. However, once nebular 
emission is taken into account (here with /esc=0) the UV spectrum becomes signif- 
icantly flatter and the trend with metallicity even inverses! In short, the UV slope 
cannot be a good metallicity indicator Figure [TT] shows the rapid evolution of /3 on 
timescales of ^ 50-100 Myr. For constant SFR, )3 reaches the typical asymptotic 
value of )3 ^ — 2 to —2.4, the precise value depending on the waveleng th base used 
to defin e j3 . Shifts of ^ 0. 1-0.3 in j8 can be typical, as illustrated e.g. by lRaiter et al 



(l20I0bh . "or PopIII we see that the nebular contribution vanishes after >20 Myr 
in a burst; between ages of ~ 10-100 Myr an integrated population can reach a 
very blue spectrum with j3 ^ —2.5-3. Otherwise, and for more extended periods 
of star formation it is difficult, if not impossible, to obtain a UV spectrum steeper 
than j3 ~ —2.5, except if a significant fraction of Lyman continuum photons are 
"leaking" (i.e. /esc > 0). Overall, the two figures show in particular the following: 
First, in principle metallicity cannot be inferred from the observed UV slope (cf. 
Schaerer & Pelidl l2005). Second, very steep/blue UV spectra (j3 <C —2.5) are not 
predicted for very metal-poor, primordial, populations, except if nebular emission is 
"suppressed". 

Some authors have indicated the possibility of unusually blue UV slopes in z ^ 7 
galaxies ( Bouwens et al. 2010l). However, the signific ance of these findings is low 



dSchaerer & de Barros 



recent measurements 



20 lot iFinkelstein et alj |2010b. and independent and more 



Castellano et al.l |20I U iBouwens et al.l 1201 It iDunlop et al 



2OI2I) do not show indications for exceptional populations or conditions (PopIII 
and strong leakage of Lyman continuum photons). Other observational aspects re- 
lated to the UV slope of distant star-forming galaxies are discussed in the Chapter 
of Dunlop. 



3.2 Ionizing photon production 

To quantity the contribution of galaxies to cosmic reionization, and to predict their 
observable spectra it is important to know the amount of energy or photons emitted 
in the Lyman continuum, i.e. at A < 912 A or at energies £ > 13.6 eV. This quantity 
is straightforwardly predicted by evolutionary synthesis models, and is generally 
expressed as a photon flux (e.g. in units of photons per second) normalized for 
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example per unit stellar mass, per unit SFR, or so. Another, observationally relevant 
way to express the Lyman continuum flux is by normalizing it to the UV output of 
the same stellar population, since the latter is a direct observable. This prediction 
is shown, for a constant star formation rate over relatively long tim escales, as a 



function of metalUcity and for different IMFs in Fig. [12] taken from iRaiter et al 
(l2010bh . 



Fig. 12 Predicted Lyman 
continuum flux as function of 
metallicity. The relative out- 
put of hydrogen ionizing pho- 
tons to UV continuum light, 
measured at 1500 A rest- 
frame, Qh/L\so{), is shown as 
a function of metallicity for 
constant star formation over 
1 Gyr. Qh/Lisoo is given in 
Li units on the left side of 
each panel, and in Ly units 
on the right. Results for dif- 
ferent IMFs are shown using 
the color codes and labels 
summarised in Tabl e 1 of 
IRaiter et alj l l2010bl) . The 
shaded area indicates the crit- 
ical metallicity range where 
the IMF is expected to change 
from a "normal" Salpeter-like 
regime to a mor e massive IM F 
(see tex t). From IRaiter et alj 
j20I0bl) . 
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As expected Qh/L\sqq increases with decreasing metallicity, since the ionizing 
flux depends very strongly on the effective stellar temperature and hence increases 
more rapidly than the UV luminosity. The IMF dependence also behaves as ex- 
pected, with the IMFs favouring the most massive stars showing also the highest the 
Qh/Li5qo ratios, since Qh increases more rapidly with stellar mass than the UV lu- 
minosity. For constant SFR and for a fixed IMF, the increase of the relative ionizing 
power from solar metallicity to Pop III is typically a factor 2 to 3. When considering 
an IMF change from Salpeter to a massive IMF (i.e. all cases except Salpeter and 
Scalo) the increase of 2///^i500 is larger, approximately 0.6 to 1 dex between solar 
and zero metallicity. For a very young (zero age) population the Lyman continuum 
production increases typically by ^40 %. 

Obviously, the ionizing photon output also d epends on age and the assumed star 



format ion history, not illustrated here (see e.g. ISchaereij 120021 l2003t : IRaiter et al 



2010bl where other normalisations are also used). For example, for zero age stellar 



populations the ionizing photon output per UV flux, Qh /iisoo- is highe r than shown 



here, typically by up to a factor 2-4 depending on the IMF (see Fig. I of lRaiter et al 
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2010bh . For very "massive" IMFs (i.e. IMFs favoring very massive stars) the pre- 
dictions shown in Fig.[T2]for constant star formation over long timescales (up to 1 
Gyr) are of course very similar to those for much shorter, probably more realistic 
timescales for the first galaxies, since in any case the lifetimes of the bulk of the 
stars is much shorter than 1 Gyr. For general ages and star formation histories these 
quantities can be derived from the available data files of the models. 



3.3 Lya emission 



Given predictions for the Lyman continuum flux from a population, the flux in the 
different H and He recombination lines and nebular continuum emission from the 
surrounding Hll region can easily be computed, as described above (Sect. 12.51 ), using 
simple case A or B recombination theory or computing detailed photoionization 
models. The equivalent widths of these lines - a very useful measure of their strength 
- can also be computed, given the stellar and nebular continuum. 

Since intrinsically the strongest line, and conveniently located in the rest frame 
UV domain, the Lya line is well known to be central to many studies of dis 



tant/pr imeval galaxies, in particular since the early work of [Partridge & Peebles 



(Il967h and the discovery of large populations of high redshift galaxies (e.g. lHu et af. 
1998t Ouchi et al. 2008 , and many other papers on such observa tions). Before the 



era of synthesis models and based on relatively crude assumptions. Partridge & PeeblesI 
(Il967h estimated that up to ^ 10% of the bolometric luminosity of primeval galaxies 
could be emitted in the Lya line. For constant star formation the most recent update, 
using our current knowledge of stellar evolution, atmospheres and the Salpeter IMF, 
place this number at ~ 3 % for solar metallicity, as shown in Fig. [13] For the same 
IMF this fraction is expected to be significantly higher at low metallicity, especially 
below Z < 1 /50 Zq, where departures from case B (due to collisional effects because 
of the high electron temperature in the nebula) significantl y increase the ioniza tion 
and hence also recombination rates, as recently showed bv lRaiter et al.l (120 1 Obi) . As 
illustrated in Fig. [13] the Lya line can therefore carry up to ~ 20-40% (depending 
on the IMF) of the bolometric luminosity in primordial gas. As already discussed 
above for the ionizing photon flux, the precise fraction of the Lya luminosity obvi- 
ously also depends on the age and star formation history. In any case, for the most 
extreme cases shown here (with IMFs dominated by massive stars), the predictions 
in Fig. [13] do not vary much with age, since the bulk of stars in such populations 
have anyway short lifetimes. 

The maximum strength of Lya, measured by its equivalent width W^ya, is rec- 
ognized as an interesting diagnostic of young, metal-poor/metal-free stellar popu- 
lations in the first galaxies. Indeed, very high equivalent widths - well beyond the 
maximum Wi^ya of 200-250 A predicted for solar metallicities - are expected 
for such populations. As such, or together with other unique spectral features of 
He II discussed below, t his line is of great interes t to search for "unusual" stel- 
lar populations (see e.g. iMalhotra & Rhoadsll2002h . While various initial predic- 
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Fig. 13 Fraction of the 
Lya luminosity to the to- 
tal bolometric luminosity, 
L(Lyo;)/Lboi for SFR=const 
as a function of metallicity 
and IMF. Results for differ- 
ent IMFs ai'e shown using 
the same colour codes as in 
Same color codes as in Figs. 
[T0l[T2] Thin lines show the 
results using standard case 
B recombination; thick lines 
the recent results accounting 
for departures from case B 
at very low metallicity. Note 
the strong increase of the 
predicted L(Lyo;)/Lboi val- 
ues from solar to very low 
metallic ity. From iRaiter et al.l 
<2010bh . 



SFR=const @ 1 Gyr 



Fig. 14 Predicted Lya equiv- 
alent width as a function of 
metallicity for very young 
(< 1-2 Myr). Thin lines 
show standard case B predic- 
tions, thick lines the predicted 
W(Lya) accounting to first 
order for departure from case 
B, leading to an increase by 
up to a factor ~ 1.5-2.5 at 
low metallicities (Z^Zcoii). 
Same color code s as in Figs. 
[MMS] From iRaiter etS] 
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tions of the Lya equivalent widths have yi elded somewhat different results (see e.g. 
Tumlinson et al .1200 it iBromm et al ]l200lh . it is now clear that non-LTE atmosphere 



models, the treatment of nebular emission, and departures from case B, are essential 
ingredients to properly predict the stre ngth and equivalent widths of Lya emiss ion 
from very metal-poor populations (see ISchaereill2002ll2003l:lRaiter et al.llioiObI) . 

With these ingredients, the maximum W]^ya predicted as a function of metallicity 
and for different IMFs reaches ^ 2000 A in the rest-frame (for a Salpeter IMF), or 
higher for IMFs favoring more massive stars, as illustrated in Fig. [14] In principle. 
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observations of such large equivalent widths of Lya, if attributable to photoioniza- 
tion from stars (as opposed to non-thermal sources), should be strong sign-posts of 
extreme conditions (metallicity and/or IMFs) expected in the first galaxies. Possible 
complications for the application of such a diagnostic include the effect of the IGM 
and dust/radiation transfer inside galaxies, which can significantly reduce the Lya 
flux (and more so than the nearby UV continuum), and hence reduce H^Lya- Further- 
more the Lya line emission varies rapidly with time (on timescales of ^ 10^ Myr). 



Neufeldl (1199 Ih has suggested that radiation transfer in a clumpy ISM could increase 



V^Lya, which, if applicable, could complicate the interpretation of high Wiya objects. 
In any case it is clear that galaxies with unusually high Lya equivalent widths (> 
200-250 A) are interesting candidates worth examine further in searches for very 
metal-poor and PopIII galaxies. 



3.4 Hardness of the ionizing spectrum and Helium line emission 
from the first galaxies 



As already mentioned above, the first stellar populations are expected to contain 
unusually hot massive stars whose ionizing spectra will be harder than that of "nor- 
mal" present-day massive stars. As a consequence these stars will emit in particular 
more ionizing photons above 54 eV (A < 228 A), the energy required to fully ionize 
Helium in the ISM surrounding these stars. This will give rise to recombination lines 
of He II observable in the UV and optical domain (e.g. He II A 1640, He II A4686), 
which should therefore be a fairly unique signature of these energetic, metal-poor 
stars. 



The use of He II lines to ident i fy PopIII stars/galaxie s was discu ssed by Tumlins on & Shull 
hOOdh : iTuniUnson et al.1 ( l200lb : Isromm et aLr ( l200lh : lOh et al.1 fcoOD; Schaerei 
(2002, 20031) . who prese nt predictions for th e strength of these lines. Schaerei 
(120031) and more recently iRaiter et al. I (l2010bl) lave discussed the transition from 
zero metallicity to present-day conditions. The predicted hardness of the ionizing 
flux of starbursts at different metallicities and the behavior of the He II A 1640 line 
as an example, are shown in Figs. [TSl and [T6l The main result is indeed a strong in- 
crease of the hardness of the ionizing flux, typically by 2 orders of magnitudes from 
metallicities Z ^ 1/50 Zq to PopIII (zero metallicity). The (maximum) equivalent 
width of He II A 1640 (again taking the "dilution" by both the stellar and nebular 
continuum into account) is generally very low. Values above > 5 A (rest frame) are 
only expected for young, very metal-poor stellar populations. One should, however, 
remember that these predictions are based on simple assumptions (case B, ionization 
bo unded regions), whic h may not always apply. 

Johnson et al. d 009l) ve computed the emission in H and He II recombination 



lines from hydrodynamic simulations taking into account the (time-dependent) leak- 
age of Lyman continuum and He+ ionizing photons. They find a smaller leakage for 
the higher energy photons, as expected, since for stellar sources these less abundant 
photons are absorbed closer to the source than H ionizing photons. They therefore 
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Fig. 15 Hardness 
e(He+)/e(H) of the He+ 
ionising flux for constant star 
formation as a function of 
metallicity (in mass fraction) 
for three different power-law 
IMFs. At metallicities above 
Z > 4. 10^"* the predictions 
from our models (crosses), 
as well as those of Leitherer 
et al. (1999, open circles), 
and Smith et al. (2002, open 
triangles) are plotted. The 
shaded area and the upper 
limit (at higher Z) indicates 
the range of the empirical 
hardness estimated from Hll 
region observations (see dis- 
cussi on in ^S chaeret 2003.) . 
From lSchaered <2003h . 
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Fig. 16 Predicted He II 
A 1640 equivalent width as 
a function of metallicity for 
very young (< 1-2 Myr) 
bursts. Same color codes as 
in Figs. [TO] [H] [13] Note that 
photoionization models pre- 
dict generally fainter He II 
A 1640 emission, hence lower 
equivalent widths, except for 
high ISM densities . From 
lRaiteretalJj2010hh . 
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conclude that the He II equivalent width should indeed be a fairly robust indicator 
for PopIII. 

As for Lya, the prediction of the intrinsic emission in Helium lines may also be 
more complicated than expected from simple "photon-counting" assumed for case 
B (and implemented in most synthesis models). The basic reason for this is that 
both H and He in the ISM compete for ionizing photons, which leads to a lower 
He+ ioniz ation rate in regions of low i onization parameter This effect, already dis- 
cussed by lStasihska & Tylenda for planetary nebulae, implies that the He II 
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lin e flux may be lower than predicted by evolutionary synthesis models, as shown 
bv lRaiter et al.l (l2010bl) . Detailed photoionization models are necessary to properly 
account for this effect. 




1620 1630 1640 1650 1660 1620 1630 1640 1650 1660 
Wavelength (A) Wavelength (A) 

Fig. 17 Left: He II A 1640 emission in the peculiar ; ~ 2 galaxy BX418 (black), compared to 
the broad and far weaker emission in the composite spectrum of 966 z ~ 2 galaxies (red). Right: 
The line is well fit by a superposition of two Gaussian components (solid red line). The broad 
component (long-dashed dark blue line) has FWHM ~ 1000 km s^', while the narrow component 
(short-dashed green line) is unresolved. The broad emission is attributed to W-R stellar winds, 
and the narrow component to nebular He II emission. The line has a total equivalent width of 
2.7 A, considera bly smaller than expectations for veiy metal-poor stellar populations. Figure from 
lErbetalJi20Tol) . 



Finally - as always - it is also useful and important to examine what is known 
empirically about He II emission. As well known, some stars (mostly the so-called 
Wolf-Rayet stars, evolved massive stars) show emission in He II lines, and these 
lines are visible in the i ntegrated spectra of some galaxies (sometimes eluded to a s 
Wolf-Rayet galaxies, cf. lContillT99lllSchaerer etal]ll999tlBrinchmann et aDl2008h . 
The He II A 1640 line is also seen in some spectra (individ ual or stacked) of high 
redshift galaxies (e.g. Shaplev et al. 2003 : Erb et al. 2010l) . as shown in Fig. (TT] 
The observed equivalent widths remain, however, small compared to the values ex- 
pected at very low metallicities. Due to the strong winds in the atmospheres of 
Wolf- Rayet stars these lines are broad, extending typically up to several thousand 



km/s jConti et alj|2008h . In spectra of sufficient signal-to-noise and spectral reso- 



lution it should therefore be possible to separate any stellar emission from nebu- 
lar He II emission (cf. Fig. [TtI i. Nebular He II emission is observed in some low- 
metallicity Hll regio ns and starburst galaxies in the nearby/low-z Universe (see e.g. 



Schaerer et al. I l 19991 for a compilation and references therein). The hardness of the 



ionizing flux inferred from the relative He II A4686/Hj3 intensities is shown by the 
shaded region in Fig. [151 it is lower by ^ 1-1.5 orders of magnitude than what is 
expected for PopIII dominated objects! However, the origin of nebular He II emis- 
sion in nearby objects remains difficult to unde rstand in many cases and several 
sources/mechanisms may contribute to it (see e.g. ISchaerei< 19961: llzotov et al I2OOII 



Thuan et al.ll2004l:lKehrig et alj|201 lUShirazi & Brinchmannll2012h . 
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Despite these open questions, the prediction of very hard spectra from primordial 
stars and stellar populations remains quite solid, with the main uncertainty probably 
being the IMF of these stars. Signatures from He II emission should therefore be a 
crucial tool for observational searches of the first stellar generations. 



4 The main observables and how to distinguish Population III? 

From the predictions from evolutionary synthesis models discussed in the previous 
Section it is already apparent how PopIII or very metal-poor populations can be 
distinguished from those of "normal" metallicities, mostly using H and He recom- 
bination lines. Other diagnostics can be derived from integrated colors, at least to 
some extent. Obviously, direct measurements of metallicity are a third way to tackle 
this question observationally. We here briefly summarize/describe these various di- 
agnostics. 



4.1 Hydrogen and Helium lines 

The predicted behavior of recombination lines of H and He+ with metallicity 
has already been discussed above. In short, unusually strong Lya emission (say 
EW{Lya)':$> 250 A) and/or strong He II emission are among the best indicators 
expected for very metal-poor stellar populations, due to their unusually high tem- 
peratures plus possibly a more "massive" IMF (i.e. favoring more massive stars). 



4.2 Metal lines 



The relative line intensities of emission line s from all elements can be predicted 
from photoionization models, such as Cloudy (iFerland et al .Il998l) . Figure[T8]shows 
the predicted spectra and strengths of several of the strongest lines from photoioniza- 
tion models as a function of metallicity. The left panel shows the overall SED includ- 
in g H and He line s for nebulae computed with blackbody spectra from calculations 
of lPanagial ( l2005h : the right panel the relative intensity of selected metal lines in the 
(rest fra me) UV and optical domain computed using realistic SEDs from synthesis 
models (Ilnouell201 ll) . As well known, the intensity of the forbidden oxygen lines 
peaks at low metallicities (Z/Zq ~ 1/50), and it is found to decrease monotonously 
to lower metallicity. Measuring e.g. a line ratio of [O III] A5007/Hj3< 0.1 would 
indicate metallicities below < 10^^ Z© (or alternatively Z/Zj >4), which could be 
feasible with N IRS pec on t he Jam es Webb Space Telescope (JWST), according to 



Panagial (l2005h and llnoud (1201 Ih . For more information on the feasibility of such 



observations see also the Chapter by Stiavelli. 




Fig. 18 Left: Predicted spectra (lines plus nebular and stellar continuum, solid) fr om photoioniza - 
tion models at different metallicities. He II and few other lines are identified. From lPanagial i2005h . 
Right: Predicted emissivities relative to H/3 of salected strong metal lines as a function of metallic- 
ity from photoionization models (filled circles). Other empirical data are shown by open symbols 
and the dashed line. From llnou3 j201ll) . 



4.3 Peculiar colors 



Several recent papers have examined/proposed various criteria which could be used 
to distinguish PopIII from other stellar populations at high redshift, based on pho- 
tometry/colors (see Inoue 2011; Zackxisson etal, 201 la b). 

As already discussed above (Sect. lTTT l. the UV slope is difficult to use as a metal- 
licity indicator due to degeneracies between age/star-formation history and metallic- 
ity, and due to the contribution of nebular continuum emission which considerably 
"dampens" the metallicity dependence of the UV slope on metallicity. In cases of a 
large escape fraction of the Lyman continuum flux the latter difficulty is of course 
not present. Predicted magnitud es and colors f or different populations a nd variable 
es cape fractions are present ed in llnoud(l201 ll) : IZackrisson et al. I (l2011bl) . 

Zackrisson et al. ( 201 Ibl) have proposed that PopIII galaxies at z ^ 8 with low 



escape fractions can be distinguished from metal-rich objects based on two colors 
between 4.4 and 7.7 jim, which could be measured with the NIRCam and MIRJ 
instrument on the the JWST. The basis for such a distinction is the presence of 
strong H lines (Ha) and the absence of [O iii] A/l 4959,5007 in these filters at these 
redshifts. 
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Fig. 19 Predicted 
JWST/NIRCAM F277W- 
F444W colors of z = 8 galax- 
ies as a function of metallicity. 
Different symbols/colors il- 
lustrate variations of the star 
formation history. Vertical 
lines connect models with 
varying Lyman continuum 
escape fraction. The horizon- 
tal dashed line is a proposed 
criterium to separate very 
metal-poor cases from higher 
metall icity cases. From ilnoud 
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Other criteria to select very metal-poor galaxies exploit the expectation of very 
blue colors between the UV and optical domain, which can to sorne ext ent be quite 
independent of nebular emission, i.e. of the escape fraction ( Inouel20"Tl ). Indeed, in 
case of strong nebular emission, the UV continuum should be strong with a Balmer 
jump in emission, and emission lines weak in the optical (due to low metallicity). 
If in contrast the escape fraction was large, one recovers again the intrinsically blue 
stellar spectrum. In both cases one therefore expects a fairly blue UV-optical color 
At z 8, for example, the color could be measured between 2.7 and 4.4 /im with 
the JWST; the corresponding predictions are shown in Fig. [19] They show how such 
a blue color could be exploited to select very metal-poor (Z<10^^ Zp) galaxies, 
although in some conditions an overlap also exists with higher metallicities. 

PopIII or very metal-poor galaxies with strong Lya emission (cf. Sect. l3.3b could 
also be detected due to the excess caused by this line. For example at z ^ 8-9 this 
can cause an unusuall y blue J — H color, as sho wn bv lPello & Schaererl (l2003b and 
discussed in detail bv lZackrisson et al.l (1201 lal) . In any case, it is clear that nebular 
emission (both line s and continua) can si g nific antly contribute to bro ad-band fluxes 
dSchaerer & de Barros 2009. .2010, 2o7lh . and lSchaerer et al.l ( 1201 lb have recently 
shown that SED modeling techniques can recover the strength of Lya emission from 
current broad-band photometric surveys of z '-^ 3-6 galaxies. This demonstrates that 
various photometric criteria and SED fitting methods should also be able to select 
peculiar objects, such as very metal-poor galaxies and PopIII dominated objects. 
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As should be clear from the onset, evolutionary synthesis models are an important, 
fundamental tool to interpret a wide variety of extra-galactic observations, from 
the nearby Universe to the most distant, first galaxies. They are commonly used to 
render complex state-of-the-art simulations of the Universe "visible", i.e. to translate 
physical properties of simulated galaxies into observables. Finally, they are also key 
for many "prospective" studies, such as for the preparation of new missions, to guide 
observers searching for Population III objects etc. 

Conceptually simple, synthesis models basically gather what is known about star- 
formation (the IMF and star-formation history), stellar evolution and atmospheres, 
and some additional emission processes, to predict the temporal evolution of the 
spectrum of an integrated stellar population. For this reason, synthesis models can 
basically only be as good/reliable as their input physics is. Therefore regular up- 
dates are necessary, and it is important to keep problems and limitations of these 
"ingredients" and assumptions in mind. 

For the spectral modeling of first galaxies, very metal-poor galaxies, and Pop- 
ulation III objects, one of the major unknowns is certainly the stellar initial mass 
function. The evolution of massive stars especially at low metallicities, observa- 
tionally inaccessible in the nearby Universe, remains also uncertain and efforts are 
ongoing to properly account for the effects of rotation, magnetic fields and related 
processes. The star formation histories of young, distant galaxies are probably com- 
plex and fairly stochastic at least on short time scales. Some of these questions have 
been addressed above, at least schematically (Sect.|2]i. 

In any case, combining the best of our knowledge evolutionary synthesis mod- 
els represent a very useful and important tool for studies of galaxy formation and 
evolution, from the nearby Universe back to first galaxies. 

Acknowledgements I have benefitted from many interesting and stimulating discussions and col- 
laborations with colleagues during the past years. I wish to thank them here collectively for this. I 
acknowledge support from the Swiss National Science Foundation. 
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